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Abstract — Quantum feedback is assuming increasingly im- 
portant role in quantum control and quantum information 
processing. In this work we analyze the application of such 
feedback techniques in eliminating decoherence in open quan- 
tum systems. In order to apply such system theoretic methods 
we first analyze the invariance properties of quadratic forms 
which corresponds to expected value of a measurement and 
present conditions for decouplability of measurement outputs of 
such time-varying open quantum systems from environmental 
effects. 

I. Introduction 

Decoherence is the process by which quantum systems 
lose their coherence information by coupling to the environ- 
ment. The quantum system entangles to the states of the en- 
vironment and the system density matrix can be diagonalized 
in a preferred basis states for the environment, dictated by the 
model of interaction hamiltonian[7][9]. Decoherence is now 
the biggest stumbling block towards exploitation of quantum 
speedup[8] using finite quantum systems in information 
processing. Many authors have addressed the control and 
suppression of decoherence in open-quantum systems by em- 
ploying a variety of open loop and feedback strategies. Effect 
of decoherence suppression under arbitrarily fast open loop 
control was studied by Viola et al [13] [15]. Another method 
along similar lines for control of decoherence by open-loop 
multipulses was studied by Uchiyama et. al.[ll]. A very 
illustrating example of decoherence of single qubit system 
used in quantum information processing and its effective 
control using pulse method was worked out by Protopopescu 
et al[22]. Shor[16] and Calderbank[17] also came up with in- 
teresting error-correction schemes for detecting and reducing 
effects of decoherence on finite quantum registers. Recently 
many authors have also studied the application of feedback 
methods in control of decoherence[14],[18]. Technological 
advances enabling manipulation, control of quantum systems 
and recent advances in quantum measurements using weak 
coupling, non-demolition principles [23] etc, has opened up 
avenues for employing feedback based control strategies for 
quantum systems [19], [20], [18]. 

In this work we analyze the effectiveness of feedback 
method in eliminating decoherence. A wave function ap- 
proach as opposed to density matrices for the schrodinger 
equation is adopted which represents the system in an input- 
affine form and greatly enables one to exploit methodologies 
from systems theory. We first analyze what it means for a 
complex scalar function to be invariant of certain parameters. 
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The generality of the treatment adopted here makes all types 
of quantum systems amenable to the results. It is also shown 
here that analysis of invariance of quadratic forms also lead 
to Decoherence Free Subspaces (DFS) for the open quantum 
systems but from a different and general perspective. DFS 
was first shown to exist by Lidar et al[12] by analysis 
of Markovian master equation for open quantum systems 
that naturally gives rise to subspaces that are immune to 
the effects of decoherence namely dissipation and loss of 
coherence. 

II. Mathematical Preliminaries 

We explore the conditions for a scalar function represented 
by a quadratic form of a time varying quantum control sys- 
tem to be invariant of perturbation or interaction hamiltonian 
when coupled to a quantum environment. 

Let 

r 

+ ^u^{t)H,®J^{t,x)]£,{t,x) 

i=l 

be the governing Schrodinger equation for a quantum 
system interacting with the environment. 
Hs be the system's Hilbert space. 
He be the environment's Hilbert space. 
Hs could be finite or infinite dimensional and He is generally 
infinite dimensional. 

^{t, x) be the wave function of the system and environment. 
Ho and He are respectively the drift Hamiltonian of the 
system and environment while Hi's are the control Hamil- 
tonian of the system. Hsb governs the interaction between 
the system and the environment. The above Hamiltonian are 
assumed to be time varying and dependent on the spatial 
variable. Consider a scalar function (typically the expected 
value of an observable) of the form, 

y{t,o^m^)\c{t,xmt,x)) (1) 

where again C{t, x) is assumed to be time-varying operator 
acting on system Hilbert space. The above is the general 
form of a time dependent quantum system and we wish to 
study the invariance properties of the function y{t,£^) with 
respect to the system dynamics. 

Let t/(t,^) ~ f{t,x,ui,---,Ur,HsB) be a complex 
scalar map of the system as a function of the control 
functions and interaction Hamiltonian over a time interval 
io < ^ < ^1- The function is said to be invariant of the 
interaction Hamiltonian if 

f{t,X,Ui, ■ ■ ■ ,Ur,HsB) = f{t,X,Ui, ■■■ ,Ur,Q) (2) 



for all admissible control functions ui, ■ • ■ ,Ur and a given 
interaction Hamiltonian Hsb- 

Let M be the manifold contained in the Hilbert space 
Tis <S> "He on which the dynamics of the system is described. 
It could be a finite or infinite dimensional submanifold of 
Sh, the unit sphere on the collective Hilbert space. The 
quantum system is assumed to be governed by time varying 
Hamiltonian and it is known that the system evolves on 
an analytic manifold V^, which is dense in Ai and a 
submanifold of the unit sphere Sh by Nelson's theorem[10]. 
Recent analysis of controllability criteria and reachability 
properties of states as studied by Schirmer et.al [30] [31] 
provides insight into behavior of quantum control systems 
on finite dimensional manifolds in R". The controllability 
under various realistic potentials was also studied by Dong 
et.al [29] [32]. However the analysis of time-varying systems 
carried out here assumes in general that the component 
Hamiltonian operators carry explicit time dependence which 
is not under the control of an external agent. And we do 
so by introducing a time invariant system in the augmented 
state space domain Ai' = A4 ® M.. A similar scheme was 
also used by Lan et.al [26] to study controllability properties 
of such time-varying quantum systems. 

Let xi — t, the new equation governing the evolution of 
the system can be written as. 



with. 
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{Ho{xi,x) + He{xi,x))(,{t,x) 


UiHi{xi,x)^{t, x) 


HsB{xi,x)^{t,x) 



y{t,0 = m,x)\C{t,xmt,x)) 



(3) 



(4) 
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The vector fields Kn — , , -^-^ t-t ,^ 
' {Ho + H,)i{x,t) 

, . , and Ki = ( c/ +^ I corre- 

sponding to drift, control and interaction can be identified to 
contribute to the dynamical evolution. 

Lemma 2.1: Consider the quantum control system (|3} and 
suppose that the corresponding output given by equation 
is invariant under given Hsb- Then for all integers 
p > and any choice of vector fields Xi , Xp in the 
set {Kq, Ki, • • • , K,.} we have 

LkjLx, ■ • • = 0; for all t, ^ (5) 

Before proving the above Lemma it is useful to consider a 
simple extension. Consider for a fixed number vector fields 
{Xi, ■ • ■ , Xp}, with p fixed and from the previous condition, 

LkiLxi ■ ■ ■ Lxj,-iLKay{t,0 =0 (^p = -^''o) 



(6) 



Combining the above conditions we get 

LkjLxi ■ ■ ■ Lxp_iLry{t,^) = 

where t G span{Ko, Ki, ■ ■ ■ ,Kr}- By finite mathemat- 
ical induction over all the variables we can replace the 
vector fields Xi , Xp with vector fields Zi, - ■ ■ ,Zp in 
spanji^o, Ki, - ■ ■ , Kr}- Hence one can show that the previ- 
ous condition is equivalent to the requirement that 

Lk.Lz, ■ --LzMt^O = 0;/or all t,C (7) 
for all p > and any choice of vector fields of the form 



Z, - Kn 



(8) 



where U , stand for the set of admissible control functions. 
Proof Now let y be invariant under Hsb- Then for small 
ti , • • • , tfe by equation (O 

y{t, Zl-oZ':zlo- - -oZ^ne)) - y{t,~Zl^o~zl^-lo. - -oZ^CO) 

where Z\, - ■ ■ ,Zk are of the form (|8j and Zi, - ■ ■ ,Zk are 
given by, 

Zk = Zk + Ki (10) 

and Zl'', the one parameter group of flow of the vector 
field Zk- The left hand side of equation (|9} is the output 
for Hsb = while the right hand side is for an arbi- 
trary Hsb- Differentiating both sides of (|9jl with respect to 
tk,tk-i, ■ ■ ■ ,ti at respectively tk = 0, ■ ■ ■ ti ~ yields, 

Lz.Lz, ■ ■ ■ LzMt.O ^ L^^L^^ - ■ - L^^yit,0 dD 
for all A: > 0. Now for k = 1 the above equation yields, 

LzMt,O^L^^y{t,0 

Since L^^^y ~ Lz^y + Lxiy, and using the above equation 
we can conclude Lxjy = 0, which is same as the equation 
for p = 0. A gain in general by induction we obtain. 



Lz.L 



■LzMt^O ^ Lz^Lz^---Lz^y{t,C) 



and using eauationfTot this yields, 

LkjLzi ■ ■■Lzt,y{t,C) = O;for all x 

for all Zi of the form (|8}. The sufficient condition for output 
invariance however requires a stronger condition of analyt- 
icity of the system. Lemma 12.11 implies that the necessary 
conditions for output invariance are. 



LK,y{t,0^^ 
LkiLk,„ ■ ■ ■ Lk: y{t, = 



(12) 



for < jQ, • • • , < r and n > 0, where A'o, • • • , Kr 
are the vector fields of the augmented system and Kj, the 
interaction vector field. The previous condition can also be 
restated thus. 



LkiLxi ■ ■ ■ Lxp_iLK^y{t,^) = {Xp = Kr) 



LKjLr,^ ■ ■■Lr^^y{t,Cj = 



(13) 



where Ti^, - ■ ■ ,Ti^ e spanjifo, ■ • • , Kr}- The above restate- 
ment might be helpful in simplifying calculations for Lie 
derivatives. 

Lemma 2.2: Suppose the system Q is analytic, then y is 
invariant under given Hsb if and only if (|5} is satisfied. 
Proof Consider a sequence of arbitrary control functions in 
U. Let 

u{t) = {u\, ■ ■ ■ ,ul),t G [tQ,ti),to = 0; 
= (ui, - •• ,ul),t e [ti,ti+t2) 



and two time instances s, t satisfying < s < i < + - • • ip. 
We can then write, 

S = <1 H h ffe_l + {tk - Tk) 

t = ti + --- + tk-i +tk + --- + t;-i + Ti 

for some index variables k. I such that < fc < Z < p 
and some r^, r/ such that < Tk < tk and < t; < i/. Let 

Xl 



'tp{t) ^ ^ ^"^^^ J be the state of the system in the augmented 

manifold and let 70 (t, s, ip{s)) and 7/(i, s, ip{s)) be the state 
map of the quantum control system in the absence and 
presence of Kj respectively, where ip{s) is the initial state at 
time s. Define a smooth function on the augmented manifold 
f{ip) — y{'yo{t, s, Making use of following relation. 



y{t, ui, U2, ■ • ■ , "r, Hsb) - y{t, ui, M2, • ■ • , Ur, 0) 



f{tp{s))*.Ki{tp{s))\.^(s)=^t(sflMo))'^^ 



(14) 



Without loss of generality, considering a piecewise constant 
control set the term inside the integral can be written as, 

/(^(s))=2/(7o(t,s,V(s))) 

= y{Zj^ o Zt-l o . . . o Z'^^l Zl" Ws))) (15) 

where Z^'s are of the form (|8}. Since the system was assumed 
to be analytic we can write, 

y{Zj^oZt-,^o...oZ'^^lZl-{^^{s))) 

7^1-1 



i 



i=0 



for some small Ti,ti^i, ■ ■ ■ such that the summation con- 
verges. The remaining terms can be expanded in the same 
way for any i, 

L'zMzti-----Zl''ms))) 

= E ^L%^-L^,^y{zl-i o . . . o Z-(^(^)) (17) 
j=o ■'■ 

Since equation MAI is zero for any t> s>Q and any given 
sequence of control functions it follows that the individual 
terms in the summation vanish yielding condition M2\ and 
hence as a consequence condition (|3 has to hold. 



III. INVARIANCE FOR THE QUANTUM SYSTEM 

Calculation of Lie derivatives The Lie derivatives in 
the above cases can be calculated for the special case 
when Ti, • • • ,Tr G {Kq^ Ki ■ ■ ■ , Kr, Kj}. For instructional 
purposes we present here two ways for calculating Lie 
derivatives of the output with respect to the vector fields 
of augmented system , 



LKiy{t) 



9(;|C(xi)|g) 3(g|C(xi)|;) 
dxi 

/ 9{i\C{x^m \ 

9e / 



■ Ki 



K 



(18) 



where Kj = ( (Cl^ss ) is the co-vector field corre- 
sponding to the vector field Ki, Hsb skew hermitian, x\, 
^* are conjugate variables and assumed to be independent of 
Xl and ^ in calculations. Therefore d{^\C{xi)\^) / dx\ — 0, 
Hence 



, Hsb\C) 



cm) 

= uc,HsBm 

Now consider, 

1 

{H^ + H,)\S) 



LK,y{t) = ( i^ic'm) mit) 
+ ( -{C\HsB 



ovit) = ( (mm) {i\c{t) 
+ ( 1 -imo 



H, 



(e|C+[C, iHo + H,)]\0 




cm) 



The variable xi is replaced with t as it was only a dummy 
variable used for calculations. Another approach follows 
directly from the geometrical interpretation of Lie derivatives 
of scalar functions, 

L^,^y(f) = lim-^(e|C(xi)|e) 
s^t ds 

with only the vector field Kt turned on for i = 
{0,1,--- ,r,/} (i.e) 

d_ / Xl 
dt \ iit,x) 
From straight forward calculations one obtains, 

d 



= K,, 



LK^yit) 

= {i\c{xi)\i) + {i\cm) + {i\c{xi)\^) 

= {i\C{t) + [C{t), (ffo + H,){t,x)m (19) 

and similarly LK,y{t) = {^\[C{t), H,{t, x)]\^) for i = 
{1, • • ■ , r, /} and Hj — Hsb- Following the above trend 
for a few Lie derivatives with respect to the vector fields 
Ko,Ki,--- ,Kr,Ki, 

LK,y^UC,HsB]\O^0 
LK,LK^y ^ mC,H,iHsB]\0 ^ 



Lk,Lk,Lk,v =mC ,H,],Hsb] 



LkiLkqLkqV = 

= LkA^\[C,Ho] + [[C,Ho],Ho] + ^C{t) 

+ |[C(t),77o(t)]|0=0 
i.e mC,H^],HsB] + [[[C,Ho],Ho\,Hsb] 
+ [C{t),HsB] + [j^[C{t),H^{t%HsBm = 

Lk.LkoLk.v = H,] + [[C, i/,], ^o]|0 = 

i.e (CI [C, H,] , i/ss] + [[[C, H,] , Ho] , i/ss] 10=0 

We are now ready to state the condition for output invariance 
of non-demolition measurements with respect to perturbation 
or interaction Hamihonian. 
Theorem 3.1: Let 



Ci = span{a(i-^^C(i)|j = 0, 1, . . . ; i = 1, . . . , r} 

Ci = |(^arfHo + ^y(?i;i = o,i,---| 

6*2 = span{ad^jj,Ci{t)\j = 0, 1, . . . ; i = 1, . . . , r} 
C2 = l (a^ffo + ^) C2;j-0,1,--- I 



C„ = span{ad^jj.Cn-i{t)\j = 0, 1, . . . ; i = 1, . . . , r} 

if , 

C„ = < I adHn + 



Similarly for = 1 we have 



"'d'Ho + tt: ) CnJ = 0, 1, 



Define a distribution of quantum operators, C{t) = 
A{Ci(i), C2(t), • • • , Cn{t), ■■■}. The output equation (□ of 
the quantum system is decoupled from the environmental 
interactions if and only if, 

[C{t),HsB{t)] = (20) 
Proof The proof follows by noting the equivalence of equa- 
tion il2\ with the above condition. Consider the following 
term Lxi^ • • ■ LKi^y{x) for any fc > 1, and io, • • • ,«fe G 
{0, • • • r}. From the calculations above it is the expected 
value of an operator of Lie brackets of Hi^ , Hi^ , • • • Hi^ , C 
and their time derivatives. In particular for fc = 0, and 

where (5(io) is the delta function that takes value 1 when 
Jo = and the operator Ti as defined is such that Ti e Ci. 



= (eiT2io 



and r2 G C2. Continuing so, in general we have r„ € C„. 
And by using condition Mil , we have [Hs b , Tn] = in 
general for decoupling. Since the condition is true for any 
n > and any T„ and since the vector space of bounded 
linear operators is complete we have [HsB,J2iLo'^i'^'i-] ~ 
Si^o ^i[HsB,Ti] — for ai e M. The converse is true by 
noting that any operator in the distribution C (i.e) for any 
T £ C can be decomposed into a sum of operators ^ cuTi 
for Ti e Ci and given [-ffss, I]i^o "^^^^l = ^^"^ which 
is true only when [HsB,Tn] = for any n. Hence from 
the previous equations LkiLk^^Lk^ _^ ■ ■ ■ Lxi^ = for 
«o, • • • ,'ik G {0, ■■■r}. 

IV. Examples 

Decoherence as studied by many authors[21][9][7], en- 
tangles the states of the system and the environment and 
amounts to forcible collapse of the wave function corre- 
sponding to preferred pointer basis decided by the environ- 
ment. The evolution of such a system can only be described 
at best at a statistical level. 

We present two qualitatively different examples to illus- 
trate the applicability of the above formalism in practical 
quantum control systems. 

A. Electro-optic Amplitude Modulation 

Consider a driven electromagnetic system in a single mode 
subject to decoherence. The control system describing the 
oscillator under the semiclassical approximation is 



UJjc]jCj + m(i)(at — a) 



where the system represented by mode a is coupled to a 
bath of infinite number of oscillators, Cj with corresponding 
coupling constants kj and where is the combined wave 
function of the system and bath. The control u{t) is the 
strength of the input current and let Hq = uja^'^a+J^j ^j'^j'^i 
and Hi = (a^ — a). Let the system be monitored by a non- 
demolition observable 

C{t) — aexp(icji) + a) exp(— iwi) 

with the corresponding output given by y{t) = 
{i>{t)\C{t)\il){t)). Following theorem 3.1, we 
have [C{t),Hi] = e'"* + e"^"^* = 2cos(wi) 
with vanishing higher order commutators. Hence 
Ci = {ci * C + C2 * I * cos(wt), Vci, C2 G K.} and 
since [C{t),Ho\ + dC/dt = we have Ci = Ci and 
the sequence converges to Ci which in general need not 



elements of the set Ci are not all zero, condition i20\ is 
not fulfilled and the non-demolition measurement is (i) 
not invariant of the interaction hamiltonian, (m) no longer 
back action evading due to the presence of the interaction. 
The measurement of the observable C{t) would thus reveal 
information about the decoherence of the system. 

B. Decoherence free subspaces(DFS) 

The techniques developed in the previous sections can be 
applied to the problem of analyzing the decoherence free sub- 
spaces(DFS) discussed in [12]. Decoherence free subspaces 
(DFS) camouflage themselves so as to be undetected by the 
interaction hamiltonian due to degeneracy of their basis states 
with respect to Hsb- 

Decoherence of a collection of 2-level systems: For a 
collection of 2-level systems interacting with a bath of 
oscillators the corresponding hamiltonian is 



converge at all. Since the commutator of the interaction Decoherence in the presence of control: In the presence of 
hamiltonian Hsb ~ + J2j '^j'^j ^^^^ '^^e the external controls = Uicr[*\ the invariance condition is 

no longer satisfied for the operator C as [[C, cj''], Cg"'''] ^ 
and hence the coherence between the states is not preserved. 
This is because of the transitions outside DFS caused by 
the control hamiltonian. The above formalism is helpful 
in analyzing general class of information that would be 
preserved in the presence of interaction hamiltonian which in 
turn would tell us about how to store information reliably in 
a quantum register Though the procedure outlined above to 
determine C (t) could get computationally intensive even for 
modest systems it is nevertheless helpful in learning about 
any ansatz C{t). 

V. Feedback Control 

The technique of using feedback has been considered 
by a number of authors [14], [19], [18] etc. Although one 
cannot extract information from a quantum system without 
disturbing it to some extent, due to rapid advances in 
quantum control technology a good deal of work carried out 
on weak measurements [6], probabilistic state estimators[l], 
non-demolition measurements and filters[28][23][25] that 
prevent systematic back action on the system, enable us 
to extract information with minimal disturbance and can 
now be applied to practical quantum systems. In reality a 
system is coupled to a probe which in turn is immersed in 
the environmental bath in order to extract state information 
of the system. The effects of feedback and probe coupling 
are currently being investigated by the authors under this 
framework. In this section we analyze the effects of minimal 
back action feedback on the control of decoherence problem 
and derive conditions for decouplability. 

Consider the augmented system equation (|3} that describes 
a time dependent quantum system and a feedback of the 
form u = a{^) + in order to preserve the input-affine 

structure of the state equation, where a, (3 are r x 1 vector 
and r X r matrix respectively of scalar functions depending 
on state |^) of the system. 



H 



N 



N 



J = l k k j=l 

where the system is assumed to interact through the col- 
lective operator J^j'^s'^ 9k ^ describe coupling to the 
mode k. An inquiry into what information about the system 
is preserved in the presence of the interaction could be 
answered by expressing the operator C acting on the system 
Hilbert space in its general form in terms of the basis 
projection operators. 



C{t) 



N)01 



E 

i,i=0..2"-l 

and solving for condition ( I20> . For a simple N=2 system we 
have after straight forward calculations 

C = span{ ^ K) (j I . (j W - z W + J - , 

id 



WK = 0,1,2...} 

where j*^'^ etc., stands for the Z*'' letter (either or 1) of the 
binary word j. Condition ( I20t . which is [C,Hsb] = now 
translates to 



_9 / XI 
at I 



(1) 



,(l)+/2) 



.^(2))^ = 0,Vi4: = 1,2,3.. 



or nontrivially, j'^^^ = i^^^ +1^2), or that the two words 

have equal number of I's. 

The above calculations are valid for any finite N, a specific 
example for iV = 3 is C = |000)(000| + |001)(001| + 
|010)(100| + |011)(101|. Of particular interest ai-e terms 
like 1011} (101 1 and |010)(100| as the coiTesponding y{t) = 
{ijj{t)\C{t)\ip{t)) which is a function of the coherence 
between the basis states |011),|101) and |010}, |100} is 
predicted to be invariant under the interaction. It is worth 
noting that the operator C{t) acting on system Hilbert space 
here need not necessarily be hermitian and only describes 
preserved information in a loose sense. 



1 

{Ho+ He+E(^rHr)ixi,x)^{t,x) 



E «i E PtjHj{xi,x)^{t, x) 


HsB{xi,x)^{t,x) 



(21) 



where again the following vector fields can be iden- 
tified Kfl = ( , T-r rT\>/ ,\ I .Kj 

° \ {Ho + H^ + Y.a^Hi)^{x,t) 

EPvmi^^t) ) ^"'^ ^' ^ ( HsBax,t) 
As stated above the necessary and sufficient conditions 
for a scalar function y{t) of the system to be invariant of the 
interaction vector field is, 



LKrV{t) = 



(22) 



for < io, - ■ ■ ,in < r and n > 0. Translating the above 
conditions into operators for the above system we obtain 
the following conditions. In the equations below we omit 
the summation symbol and following Einstein's convention 
a summation has to be assumed where ever a pair of the 
same index appears. 



L^L^^y={^\[C,f3uHi 



[[C,H + a,HjlPuH,] 
+ [C,H,]Lji^a,\O=0 

= mC,PaHi],HsB] + [[C,Hj]Lj,a,,HsB] 
+ [[[C,Ho + ajHj],PaHilHsB] 

+ [[C, Ho],Hi]LK,Pri + [[C, H,],Hi]LK,a,Pa\0 
=0 (23) 

The first two lines of RHS of the above equality is found to 
belong to the distribution [C{t), Hsb] and the last two lines 
belong to C{t). The above calculation can be extended to 
any number of terms to encompass the result. In general one 
finds that, in the presence of feedback terms the condition 
for decouplability is relaxed to 

[Cit),HsB]cC{t) (24) 

In order to solve eq.(l23> and consequently J24> for the 
feedback parameters, it has to be noted that the first two 
lines and last two lines of eq.(l23> denote operators acting 
on different Hilbert spaces, namely the system-environment 
and just the system respectively and the two terms cannot 
be reconciled unless they vanish individually which leads us 
back to original conditions for open loop invariance. 

In other words, in order for the feedback to be an ef- 
fective tool in solving the decoherence problem, the control 
hamiltonians Hi have to act non-trivially on both the Hilbert 
spaces which would enable all the operators in ( I23> act on 
system-environment Hilbert space. 

VI. Invariant Subspace 

As stated above the fundamental conditions for invariance 
were, 

LkjLk,^ ■ ■ ■ LKi^^ y{t, = 

where < ia, ■ ■ ■ j^n < rin > 0. We now explore a larger 
class of vector fields AV containing Kj that also satisfy the 
above conditions, i.e. 



(25) 



Set of such vector fields form a vector space or a distribution 
and constitute a invariant distribution in the sense described 
by the following theorems. 

Definition The vector field Kr satisfying equations 
(125 > is said to be in the orthogonal subspace of 
the observation space spanned by the co-vector fields 
dy{t,^),dLK,g ■ ■ ■ LK,^y{t,£,), ■■ ■ for all < «o, < 
r and n > 0. Denoted by Kr G 

Lemma 6.1: The distribution is invariant with respect 
to the vector fields i^To i • ' ' i Kr under the Lie bracket opera- 
tion, (i.e) if Kr e O^, then [K^, K,] € for i = 0, ■ ■ ■ ,r 
proof: Assuming a form for the vector field Kr = 

jj the Lie bracket of [Kr,Ki] for i — I,-- - ,r 

can be computed as follows. 



[Kr,K,] = 





Hr\0 Hr 





H^\^) 







Hr\i) 



[Hr,HM) 



Now using Jacobi identity. 



= -LK\LKiy{t) - LK,LK,y{t) 
= 



Now for i = and Ka = 



1 



we have. 



[Kr^Ko] 


Ho 



(26) 









Hr\0 Hr 



1 

Ho\0 



{[Hr,Ho]-Hrm 
L[Kr,Ko]yit) = imC, [Hr,Ho]] - [C,Hrm 



(27) 



We already have, 

LKrLK„y{t,£,) 

LKoLK^y{t,(,) 



{S,\[C,Hr] + [[C,Ho],Hrm^Q 

{S,\j^[C,Hr] + [[C,Hr],Hom 



Adding the above equations and using Jacobi Identity we 
conclude that [Kr,Ko\ e O^. 

VII. Conclusion 

We analyzed the conditions for eliminating the effects of 
decoherence on quantum system whose coherence can be 
monitored in the form of a scalar output equation. The results 
hold globally on the analytic manifold. 



VIII. Future Work 

The invariant distributions possess many desirable qual- 
ities and helps in control of decoherence. We wish to 
construct an algorithm to determine the invariant distribution 
for a given quantum system and its interactions. Design and 
study of feedback and analysis of the resulting stability for 
quantum control system will help us solve the decoherence 
problem for practical quantum systems. The results can be 
extended and conditions can be derived for different types 
of measurements and information extraction schemes. 
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